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Abstract —The results are given of experimental investigations of the kinetics of nucleation in gas-saturated 
helium-oxygen solutions. The temperature dependences of the mean lifetime of solution in the metastable 
(superheated) state are determined for fixed values of pressure of 1.171 and 1.667 MPa and for several values 
of concentration. The experimental data are compared with the theory of homogeneous nucleation in a macro¬ 
scopic approximation with due regard for the dependence of the properties of the nuclei on their size. The 
dimensional effect is included within van der Waals’ gradient expansion. The problems associated with the con¬ 
struction of the kinetic theory of boiling of gas-saturated solutions are discussed, and good agreement is ascer¬ 
tained between the experimental data and the theory of nucleation when the latter includes the dimensional 
effect. 


INTRODUCTION 

A supersaturated solution of gas in liquid is an 
example of a two-component metastable system in 
which the nonequilibrium is caused by the excess of the 
concentration of the dissolved (gas) component over 
the value of concentration in the saturated state. The 
supersaturation is removed as a result of isolation of the 
gas phase in the form of numerous bubbles (“cham¬ 
pagne effect”). 

A special situation is observed for a gasified liquid 
at temperatures close to the critical point of the solvent. 
Here, a two-component system may prove to be “dou¬ 
bly” metastable. Given the temperature T and the con¬ 
centration of dissolved component c\ a solution at a 
pressure p below the saturation pressure p s is supersat¬ 
urated and characterized by the degree of metastability 
A p(T, c') = p s - p, with the pure solvent being super¬ 
heated by A77 p) = T-Tf c' = 0), where T s (c = 0) is the 
saturation temperature of pure solvent at a given pres¬ 
sure. 

Experimental investigations of limiting superheat¬ 
ing AT(p, c') = T - T s of solutions of carbon dioxide and 
nitrogen in organic liquids at a low nucleation rate (J — 
10 7 s -1 nr 3 ) were performed by Mori et al. [1], by For¬ 
est and Ward [2], and, at a higher nucleation rate (J — 
10 24 —10 26 s -1 nr 3 ), by Skripov and Pavlov [3]. Mori 
et al. [1] and Forest and Ward [2] used the method of 
floating-up droplets, and Skripov and Pavlov [3] 
employed the method of pulsed heating of a wire probe. 
The gas saturation of liquid resulted in a shift of the 
boundary of spontaneous boiling towards lower tem¬ 
peratures. The authors of [1-3] observe an adequate 
agreement between the classical theory of homoge¬ 
neous nucleation and the experimental data on the 
value of limiting superheating. 

This study is aimed at determining the temperature, 
baric, and concentration dependences of the rate of 


nucleation in gas-saturated solutions. The employed 
method of measuring the lifetime of metastable liquid 
enables one to determine the temperature of superheat¬ 
ing the solution with an accuracy an order of magnitude 
higher than that of the methods used in [1-3]. The 
experiments are performed with solutions of helium in 
liquid oxygen. Gas-saturated solutions of cryogenic 
liquids belong to the class of weak solutions. At a pres¬ 
sure of 5 MPa, the solubility of helium in oxygen does 
not exceed 2-3 mol % [4]. A significant difference of 
the parameters of the potential of molecular interaction 
between the solvent and solute leads to a considerable 
nonideality of the two-component system; in this case, 
helium behaves as a surfactant and causes a decrease in 
the excess free energy of the liquid-vapor interface [5]. 
The simplicity of the molecular structure of solutions 
enables one to use van der Waals’ theory of capillarity 
for including both the adsorption and dimensional 
effects in nucleation. 


EXPERIMENT 

The experimental setup and procedure were 
described in [6]. A gas mixture of preassigned concen¬ 
tration was condensed into a measuring cell in the form 
of a glass capillary sealed on top and communicated 
with a metal sylphon bellows. The capillary was ther¬ 
mostatically controlled at a temperature T, and the tem¬ 
perature of the sylphon bellows was maintained close 
to that of normal boiling of oxygen. The supersatura¬ 
tion of the solution in the capillary was developed by 
releasing the pressure to a preassigned value of p. The 
concentration c' of the gas component in the liquid 
solution was determined by the data on phase equilib¬ 
rium [4] using the method described in [6]. In the 
experiment, the time x from the moment of transforma¬ 
tion of the solution to the metastable state to its boiling 
was measured for fixed values of T,p, and c'. After boil- 
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ing, the pressure in the system was increased to a value 
of p b > p s , where p s is the saturation pressure. The solu¬ 
tion was “held” in this state for some time z b , after 
which the measurements were repeated. The mean time 
of waiting for boiling x = Z,x,//V was found by the 
results of IV = 20-100 measurements. The error of 
determination of temperature was 0.05 K; that of pres¬ 
sure. 0.01 MPa; of time, 0.01 s; and of concentration, 
0.02 mol %. 

The experimentally recorded time X includes the 
time of waiting for the emergence of the first critical 
bubble, x w , and the time x,, of its subsequent growth to a 
macroscopic size (x = x w + x„). In the investigated range 
of parameters of state of helium-oxygen solutions, T g — 
10 3 s. The characteristic time of experiment is at least 
0.1 s. The condition x > X,, enables one to assume x = x w 
and calculate the nucleation rate as / = (x V) -1 , where 
V = 68 mm 3 is the volume of metastable liquid. 

The experiments were performed for two values of 
p- 1.171 and 1.667 MPa and several values of concen¬ 
tration of helium in oxygen. The investigated sub¬ 
stances were 99.99% (He) and 99.92% (0 2 ) pure. The 
values of the mean lifetime of metastable solution 
being measured were in the range from 0.1 to 1000 s, 
this corresponding to the nucleation rate values of 
10 4 -10 8 s" 1 nr 3 . 

In the case of fixed external conditions (T, p, c - 
const), the emergence of the first viable nucleus in the 
metastable liquid is an event of the Poisson type [7], 
A typical histogram of the experiments (the distribu¬ 
tion of the number of boilings i falling within the inter¬ 
val of x, X + Ax) is given in Fig. 1. The experiments 
revealed the dependence of the measurement result on 
the “hold” time x h . For the values of time x h — 3 min, 
which are characteristic of experiments with one-com¬ 
ponent systems, boilings were often observed with the 
values of hold time that were “anomalously” small for 
the preassigned values of T, p, and c'. “Anomalously” 
small values of x could appear and disappear in the 
course of experiment or be retained throughout the 
measurements. This behavior of a gas-saturated system 
is analogous to that of a one-component system, when 
readily activated boiling centers are present in the latter 
[8], The law of distribution of “anomalously” low val¬ 
ues of the time of waiting for boiling is likewise close 
to Poisson’s law (Fig. 1). When the “hold” time 
between measurements increases to 30-40 min (that is, 
by an order of magnitude!), no premature boiling of liq¬ 
uid was observed. In the course of processing of exper¬ 
imental data, the “anomalously” low values of x were 
eliminated from treatment. 

Figure 2 gives, on a semilog scale, the temperature 
dependence of the nucleation rate in helium-oxygen at 
fixed values of pressure. The vertical dimension of the 
points corresponds to the statistical error of determina¬ 
tion of J. As in the case of one-component liquids, two 
segments with different behavior of the temperature 



Fig. 1. A histogram of experiments (c' = 0.14 mol %, 
p = 1.171 MPa, T = 166.66 K): Solid lines indicate N = 29, 

X = 2.66 s, and x b = 40 min; broken lines indicate N = 66, 
x = 1.34 s, and x b = 3 min; and smooth curves indicate the 
results of calculation by the exponential law. 

dependence of the nucleation rate may be identified on 
the J(T) curves. In the case of low superheating, the 
experimental isobars have characteristic bends; as the 
temperature rises, these bends are replaced by segments 
with a sharp increase in the nucleation rate. The first 
segment is associated with initiated nucleation, and the 
second one is associated with spontaneous boiling of 
liquid solution. The dissolving of 0.1 mol % of helium 
in liquid oxygen causes an approximately 10% reduc¬ 
tion of the limiting superheating A T sh (p, c) = T s h ~ T s in 
the region of spontaneous boiling of liquid. In so doing, 
the mean lifetime of solution in the region of initiated 
nucleation is reduced by a factor of three-four. Note 
that the gas saturation in this case shows up at very low 
concentrations of helium. 

High-energy particles, generated by cosmic radia¬ 
tion and radiation background, are usually treated as 
the main factor which initiates the boiling of liquid and 
leads to the bending of isobars [7], However, it is diffi¬ 
cult to explain so strong an effect of gas saturation on 
the decrease in the stability of metastable solution in 
the initiation region while remaining within the frame¬ 
work of this model alone. 

In the regions where the nucleation rate increases 
sharply, the magnitude of the derivative d In J/dT within 
the experimental error does not depend on concentra¬ 
tion. At p= 1.171 MPa, it is ~18 and, at p= 1.667 MPa, 
it reaches 21. This corresponds to the increase in the 
nucleation rate by approximately eight-nine orders of 
magnitude with a temperature rise by 1 K. 
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Fig. 2. The temperature dependence of the nucleation rate at a pressure of (a) 1.171 MPa and (b) 1.667 MPa and concentration of 
(1) c' = 0, (2) 0.08, (3) 0.14, and (4) 0.20 mol %. Dashed line indicates the results of calculation by the theory of homogeneous 
nucleation with a = (oxygen), and dotted line indicates the same but with a = c(R*) (oxygen and helium-oxygen solution). 


NUCLEATION THEORY 

The theory of boiling of gasified nonvolatile liquids 
such as solutions of gas in liquid was treated by Derya- 
gin and Prokhorov [9] and by Kuni et al. [10]. In these 
studies, the thermodynamic state of the nucleus was 
characterized by two variables, namely, the bubble vol¬ 
ume v and the pressure of the gas component in this 
bubble, p ". The general case of boiling of a metastable 
solution is treated in [11], where, according to the 
Zel’dovich-Kramers scheme [12], the fluctuation 
growth of near-critical bubbles is presented as a diffu¬ 
sion process in the field of thermodynamic forces in the 
space of the variables such as the bubble volume v and 
the partial pressure of mixture components in the bub¬ 
ble p" and p '2 . The classical expression 


components of absorptive-emissive interaction 
between the bubble and solution, namely, the diffusive 
delivery of substance to the bubble and the free-molec¬ 
ular exchange at the bubble boundary, are treated as 
independent processes occurring in parallel with each 
other, the rate of variation of the number of molecules 
of mixture components in the bubble will be defined by 
the expression 


Ni = (NuH + NiDy 1 , N 2 = N lm , (2) 


where 


Nu 


kT 


(Pi.-Pi) 


(3) 



C v £exp-^* (1) 


is the rate of molecular exchange; 


was derived for the steady-state nucleation rate. Here, 
C v is the normalization constant of the equilibrium vol¬ 
ume distribution function of nuclei; R... and W... denote 

the radius and work of formation of the critical bubble, 
respectively; c is the surface tension of the bubble; X 0 is 
the decrement of the rise of an unstable variable 
(nucleus volume); and B is a kinetic factor. The dynam¬ 
ics of bubble growth are included mainly via the decre¬ 
ment Xfl for which a cubic equation was derived in [11], 
In solving the dynamic problem, it was assumed that 
the motion of the molecules of the gas component in the 
solvent towards the bubble was diffusive; however, the 
exchange of molecules between the solution and bub¬ 
ble occurred only in a free-molecular mode. If two 


Nid = 4nRp'D g H(p'U-p") (4) 

is the rate of diffusive delivery; [3, is the condensation 
number; v ti = (tikT/nm,) 111 denotes the mean thermal 
velocity of motion of molecules of mixture components 
in the bubble; m, is the molecular mass; p' is the density 
of the solution; D g is the diffusion coefficient of the gas 
component in liquid; and H is Henry’s law constant. 

Much as this was done in [11], one can demonstrate 
that, with such treatment of interaction between the 
bubble and solution, the diffusion tensor of nuclei in the 
space of their essential variables will now be symmet¬ 
ric, and the decrement of the rise of the unstable vari- 
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able will be determined by solving the equation 
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and r) is the viscosity of the solution. 

If it is required that p 2t . —*- 0, (3 2 —► 0, and c' 0, 
which corresponds to the boiling of gasified nonvolatile 
liquid, the cubic equation (5) reduces to quadratic. The 
diffusive delivery of substance to the bubble may be 
ignored if the characteristic time of diffusion is much 
less than the other time scales associated with the vari¬ 
ation of pressure in the bubble. In this case, V —► °°, and 
Eq. (5) yields the equation in [9] for determining X () in 
a viscous nonvolatile solution in which the concentra¬ 
tion is uniform up to the separation surface of the bub¬ 
ble. If the free-molecular mode of evaporation of gas 
into the bubble is ignored, i.e., —*- 0, Eq. (5) 

changes to the equation in [10] and, simultaneously 
with Eq. (1), defines the nucleation rate in a nonvolatile 
viscous gasified liquid, when the transport of gas mol¬ 
ecules to the bubble is limited by diffusion alone. 

The calculation results demonstrate that, in all 
approximations, the value of the kinetic factor B in for¬ 
mula (1) depends very little on temperature, pressure, 
and concentration, compared with the exponential fac¬ 
tor. Neglect of the volatility of the solvent and its vis¬ 
cosity leads to an increase in B compared with its value 
calculated from the general formula (5). Now, neglect 
of the diffusion mode of delivery of substance to the 
bubble likewise causes an increase in the value of B. In 
the region of spontaneous boiling of a helium-oxygen 
solution, the difference between the values of B calcu¬ 
lated for different limiting cases, when the dynamics of 
bubble growth are limited by one or another factor, 
does not exceed one order of magnitude of B. This 
causes uncertainty in the temperature of limiting super¬ 
heating of liquid of the order of 0.1 K. A si mi lar uncer¬ 
tainty in the superheating temperature is introduced by 
the inaccuracy of the order of 3% in the work of forma¬ 
tion of a nucleus. 


WORK OF FORMATION AND SURFACE 
TENSION OF CRITICAL BUBBLE 

In Gibbs’ method of separation surfaces, we have, 
for the work of formation of a critical nucleus [7], 


W* 


16 7tQ 3 
3 (pl-p'Y 


( 6 ) 


The pressure in the nucleus p* is defined by the 
condition of its material (chemical) equilibrium with 
the surrounding metastable phase, 

p*,c*) = p',c'), (7) 


and the correlation between the surface tension a and 
the radius of the tension surface of the bubble is given 
by Laplace’s relation, 

pi ~ p = 2 g/R*. (8) 

Here, p, is the chemical potential of the /th compo¬ 
nent. 

A rigorous approach to the calculation of W* calls 
for the knowledge of a as the interface curvature func¬ 
tion and of the equation of state describing the stable 
and metastable regions of the solution. The problem on 
finding c(R*) and on calculating W* without introduc¬ 
ing the concepts of separation surfaces may be solved 
if Gibbs’ method is supplemented with van der Waals’ 
method of gradient expansion [13]. According to the 
latter, the increment of a high thermodynamic potential 
during the formation of spherical nonuniformity in a 
binary system is given by the expression 


AQ{p!,p 2 } = 47t | 


o L 


* V 7 dp jdp i 


U j = 1 


r dr , 


(9) 


where 


2 

aco = p’-p + ^[p,(p„p 2 )-p;]p, 

i = 1 

2 

= /(pnp2)-/'-£(p,-p;)^ 

i = 1 

/ is the density of free energy of a homogeneous sys¬ 
tem, and kjj is the matrix of influence coefficients. 
A single prime indicates the initial value of a parameter 
in the liquid phase, and a parameter without a prime 
indicates the running local value. 

From the conditions of stability of a two-phase two- 
component system follows 

k\\k 22 — ^12 > kn>0, k 22 > 0. (10) 


If 


k l2 = (k n k 22 ) m (11) 

is assumed for the cross coefficient k l2 [14], inequality 
(10) changes to an equality, which conforms to the 
indifferent equilibrium of the system relative to density 
gradients opposite in sign and defines the limiting per- 
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missible value of the coefficient k l2 . In view of (11), 
expression (9) takes a simpler form. 


close-to-ideal gas. This enables one to use the equation 
of state for oxygen for describing the helium properties, 


G{pi. P 2 } 



( 12 ) 


Here, p p = (k n /k 22 ) 1,2 p l + p 2 . The distribution of 
density p, (i = (3, 1, 2) in the critical bubble is found 
from the solution of a set of Euler equations for the 
functional (12), which includes one differential equa¬ 
tion 



(13) 


with the boundary conditions of pp —*- p p at r —*• oo 
and dp^/dr —*- 0 at r —► 0 and r —* oo and one alge¬ 
braic equation 

Pi —Pi = ( k n /k 22 ) (p 2 —p 2 )- (14) 


In so doing, the work of formation of a critical bub¬ 
ble is 


W * = minmaxA£2{pip 2 }. (15) 

Expression (15) provides a way toward the calcula¬ 
tion of W* that is alternative to (6) and, simultaneously 
with (6)-(8), enables one to determine the correlation 
c(R *). The solution of this problem implies the knowl¬ 
edge of the equation of state for binary solution and of 
the influence coefficients ky. The values of ky are func¬ 
tions of temperature and depend little on density 
[13, 14], One can ignore the dependence of ky on p, and 
determine the coefficient k 22 for pure solvent by the 
data on surface tension at a plane interface, 

p; 

= 2k 22 2 |(Aco) 1/2 r/p. (16) 

p.l 

In so doing, the distribution of density pp in the tran¬ 
sition layer is found from the solution of the equation 

m ?) 2 = a “ < 17) 

with known boundary conditions. 

The equation of state for the helium-oxygen sys¬ 
tem, describing both the stable and metastable regions, 
was derived within the one-liquid model of solution 
based on the unified equation of state for oxygen in 
[15], Because the range of existence of the liquid phase 
for oxygen lies much higher than the critical tempera¬ 
ture of helium, helium in this case may be treated as a 


p/pRT = 1 + :£/vpT \ (18) 

i.j 

with the correction factor C, introduced into the latter 
equation. 

Here, p = p/p c2 ; T = T/T c2 , p c2 and T cl denote the 
density and temperature at the critical point of oxygen, 
respectively; and R is the universal gas constant; the 
values of the coefficients by are given in [ 15 ], In pure 
oxygen, £ = 1; at £ = 0.431, the equation of state (18) 
with the coefficients by for oxygen provides the most 
optimal description of the p, p, T properties of helium 
in the range of reduced density 0 < p < 1 and tempera¬ 
ture 0.6 < T <1. 

In describing the helium-oxygen solution, three 
individual parameters of Eq. (18) T c , p c = l/v c , and £ are 
treated as functions of concentration, 

T c (c) = T cl c 2 + 2a T T cn c(l-c) + T c2 (l-c)\ (19) 
v c (c) = v cl c 2 + 2a v v cU c{\-c) + v c2 (l-c) 2 , (20) 


C(c) = Cffi 2 + 2a c C I2 c(l-c) + (l-c) 2 , 

(21) 

where the cross terms T cl2 , v cl2 , and £ 12 are given by the 
relations 

rwi / rri r 1 1 \ 1 / 2 

1 cl2 = y 1 cl 1 cl) ’ 

(22) 

, 1/3 , 1/3,3 

v cl2 = ( v c 1 + v c2 ) /8, 

(23) 

c4 

vjp 

+ 

II 

(N 

vjp 

(24) 

The free parameters of Eqs. (19)—(21) a T , 

a v , and 


oq may be determined by the saturated vapor pressure 
and by the composition of coexisting phases of the 
solution. However, the available literature data on p s , 

c' , and c" [4] proved to be insufficient for reliable 
reproduction of the phase diagram in the temperature 
region adjoining the critical point of the solvent; there¬ 
fore, the data about the surface tension of solution [5] 
was used as additional information. This enabled us 
both to improve the reliability of description of the 
p, p, T, c parameters and to determine the influence 
coefficient of pure helium k n . The calculations were 
performed by Eqs. (18)—(24) in [16, 17]. The data on 
the surface tension of pure oxygen were borrowed from 
[16]. The results of determining the free parameters of 
the unified equation of state for helium-oxygen solu¬ 
tion are given in [17]. 

The distribution of the density of mixture compo¬ 
nents in a critical bubble, obtained by numerical inte¬ 
gration of the set of equations (13) and (14), is given in 
Fig. 3. Unlike solutions with complete solubility of 
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components (of the Ar-Kr type), the density of the first 
(gas) component in the bubble is higher than in the sol¬ 
vent. The presence of a weak maximum on the density 
profile on the gas phase side points to excess adsoiption 
of the gas component in the surface layer of the 
nucleus. 

The tension surface radius R * and the surface ten¬ 
sion a of the critical bubble are defined by the set of 
equations (6)-(8) and (15). For the radii R e ,• of equirno- 
lecular separation surfaces of the components, we have 


= & 


1/3 


. A P/ = p;-p;', 


(25) 


h = 47tJ| Pf -p; 


'I rdr. 


The behavior of the Tolman parameters of the com¬ 
ponents of binary solution (8, = R e , -/?* , = 1, 2) is sim¬ 
ilar to that of the Tolman parameter 8 of pure sub¬ 
stances: they increase monotonically with the curvature 
of the interface. In the region of large curvature radii (in 
the region of weak metastability), the relationship 
8,(1//?*) is linear. In the case of small radii of the ten¬ 
sion surface (region of strong metastability), the Tol¬ 
man parameters diverge as T?* 3 . As distinct from 8, and 
8 2 , the quantity 8 1(2) = 8, — S 2 = /?,. i - R e , 2 demonstrates 
a weak dependence on both the curvature of the separa¬ 
tion surface and the concentration. 

Figure 4 gives the surface tension of critical bubbles 
of a helium-oxygen solution in the region of spontane¬ 
ous boiling as a function of curvature of the tension sur¬ 
face. Minor additions of helium to liquid oxygen do not 
change the behavior of the dimensional dependence of 
surface tension of the solvent, while decreasing this 
tension monotonically. The inclusion of the depen¬ 
dence <r(/?*) in the work of critical nucleation causes a 
decrease in the height of the activation barrier and facil¬ 
itates the boiling of the solution. 


COMPARISON OF THEORY 
AND EXPERIMENT 

Solution (1) is valid after some incubation period of 
time t, following which the nuclei, irrespective of their 
initial distribution, are generated with a regular fre¬ 
quency J. The time t, may be estimated as the time of 
diffusive climbing of nuclei of the new phase over the 
critical region of the potential barrier by the unstable 
variable t, = t v - 1/|A 0 |. As was demonstrated in [18], for 
helium-oxygen solution in the region of parameters of 
state in which the experiments were performed, the 
time t l — 10 x s, this justifying the neglect of nonstation- 
arity. 

The exponential term is predominant in the expres¬ 
sion for nucleation rate (1). The value of kinetic factor 
B depends little on temperature, pressure, and concen¬ 


p, kmol m 3 



Fig. 3. The distribution of the density of mixture compo¬ 
nents in a critical bubble (T = 140 K, c' = 0.1 mol %, p = 
1.0 MPa). 


G, mN m 1 



Fig. 4. The surface tension of critical bubbles in a He-0 2 
mixture at T = 140 K and concentration of c’ = ( 1 ) 0, (2) 0.1, 
and (3) 0.2 mol %. 


tration. Because of the very strong dependence of the 
nucleation rate on the exponent, the differences in the 
estimates of the coefficient B are of little importance. 
The indeterminacy of the value of the kinetic factor B 
of one-two orders of magnitude leads to indeterminacy 
of 0.1-0.2 in the temperature of attainable superheating 
of liquid, which does not exceed the error of its experi¬ 
mental determination. 

Another important aspect associated with the use of 
formula (1) is the value of the normalization constant of 


HIGH TEMPERATURE Vol. 38 No. 6 2000 






858 


BAIDAKOV, KAVERIN 


Table 1 . The temperature of attainable superheating of helium-oxygen solutions, J - 10 7 s 1 m 3 


c\ mol % 

T s , K 

p s , MPa 

T s fr K 

experiment 

theory* 

theory** 



p= 1.171 MPa 



0 

122.40 

2.570 

138.17 

138.43 

137.96 

0.08 

118.96 

2.640 

137.79 

138.07 

137.58 

0.14 

114.55 

2.686 

137.43 

137.80 

137.28 

0.20 

- 

2.732 

137.06 

137.52 

136.99 



p = 1.6< 

57 MPa 



0 

129.05 

2.786 

140.00 

140.28 

139.91 

0.08 

127.12 

2.860 

139.72 

139.91 

139.53 

0.14 

125.34 

2.910 

139.46 

139.63 

139.23 

0.20 

123.04 

2.963 

139.23 

139.35 

138.94 


* Calculations for a = by the data of [5], 
** Calculations for a = o(R*), our results. 


the equilibrium volume distribution function of nuclei, 
C v . It is usually assumed that a nucleus may form on 

any molecule of the system; then, C v = p' = p+ pi, 
where pand p, denote the number density of solvent 
and solute in the liquid phase, respectively. A more rig¬ 
orous approach to the determination of B is given by 
Deryagin et al. [19], according to which C v = v*p' 2 . 
The transition from the normalization constant C v = p' 
to C v = v*p' 2 results in a decrease in the temperature of 
attainable superheating of helium-oxygen solution by 
approximately 0.2-0.3 K. In what follows, in compar¬ 
ing theory and experiment, we use the normalization 
constant for the nucleus distribution function suggested 
by Deryagin et al. [19], 

Table 1 gives data on the temperature of attainable 
superheating of helium-oxygen solutions, which corre¬ 
spond to the nucleation rate J = 10 7 s" 1 nr 3 . The exper¬ 
imental results are compared with the results of calcu¬ 
lation of T n by the theory of homogeneous nucleation in 
a macroscopic approximation (a = <7,J and with due 
regard for the dependence of the bubble surface tension 
on the curvature of its separation surface (a = 0(7?*)). 

The kinetic factor B was determined by the decrement 
found from the solution of Eq. (5). Figure 2 gives the 
temperature dependence of J in view of <7(7?*), calcu¬ 
lated by the theory of homogeneous nucleation and, for 
pure oxygen, in a macroscopic approximation as well. 
The data of Tables 1 and 2 demonstrate good agreement 
between experiment and theory, when the latter 
includes the dimensional effect. The divergence 
between the theoretical and experimental values of the 
temperature of attainable superheating is within the 
limits of the experimental error and the accuracy of the 
thermal parameters appearing in (1). This agreement is 
an indirect vindication of the theoretical assumption of 


the equilibrium composition of vapor-gas mixture in 
the critical nucleus. 

The typical parameters of critical bubbles under 
experimental conditions are given in Table 2. With a 
temperature T = 138 K and pressure p = 1.171 MPa, the 
radius of critical bubble in superheated pure oxygen 
f?* = 3.83 nm and the number of molecules in the bub¬ 
ble — 485. The dissolution of 0.2 mol % helium in 
the bubble increases the critical bubble size to f?* = 
4.26 nm. In so doing, the concentration of helium in the 
bubble c* = 6.8 mol %, and the number of molecules 
77 * — 640, and the bubble is in fact free of helium 
atoms: 77”* — 50. Before the formation of the bubble, 
its volume was occupied by 77 \ = 5674 molecules of 
oxygen and n\ =13 atoms of helium. In order to attain 
an equilibrium composition, the lacking =37 atoms 
may be delivered to the growing bubble in the process 
of chaotic thermal motion. Therefore, the process of 
nucleation in a gas-saturated solution at a temperature 
close to the temperature of attainable superheating of 
solvent and to its critical point differs little from the 
process of nucleation in a superheated pure liquid. 
A decrease in the temperature of superheating of liquid 
as the gas component is dissolved in the liquid is prima¬ 
rily associated with the shift of the line of equilibrium 
coexistence of phases and, to a smaller extent, with the 
valuation in surface tension. 

The situation changes drastically if the nucleation 
occurs at low temperatures. For example, at a tempera¬ 
ture T = 100 K and pressure p = 1.171 MPa, when the 
concentration of helium in oxygen is c' = 0.03 mol % 
and the solution is supersaturated so that the Gibbs 
number G = W*/kT - 80, i.e., close to that which we 
had in experiments in the boiling of helium-oxygen 
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Table 2. The typical parameters of critical bubbles of pure 
oxygen and of helium-oxygen solutions at a pressure p = 
1.171 MPa in the vicinity of the boundary of spontaneous 
boiling 


Parameter 

c', mol % 

0 

0.1 

0.2 

T, K 

138 

137 

136 

p s , MPa 

2.550 

2.583 

2.626 

W*/kT, Eq. (6) 

89.78 

107.09 

124.80 

W*/kT, Eq. (15) 

74.78 

91.39 

108.28 

log/, Eqs. (1), (5), and (6) 

1.07 

-6.89 

-14.83 

log/, Eqs. (1), (5), and (15) 

7.55 

-0.07 

-7.65 

a TO , mN/m 

2.47 

2.63 

2.79 

j, nm 

- 

0.289 

0.276 

8 00j 2 , nm 

0.0363 

0.0367 

0.0373 

R*, nm 

3.828 

4.066 

4.262 

o(R*), mN/m 

2.32 

2.50 

2.67 

»* 

485 

566 

640 

n 


22 

50 

i* 




5 b nm 

- 

0.429 

0.403 

S 2 , nm 

0.209 

0.188 

0.172 

c* , mol % 

0 

3.25 

6.80 


solutions at high temperatures, the critical bubble 
radius will be R* — 1.7 nm and the number of helium 

molecules in the bubble n = 146, with the number of 
oxygen molecules — 3. Prior to nucleation, 

approximately 12 atoms of helium were present in a 
volume of critical size. For an equilibrium to set in, the 
gas must diffuse to a distance of approximately 
2.5 times the bubble radius. In this case, the character¬ 
istic time of diffusion is x D = (2.5R*) 2 D^ 1 —2x 1CL 8 s, 
which is comparable with the time of relaxation to 
steady-state size distribution of nuclei. In addition, a 
stronger dependence of the surface tension of critical 
bubbles is possible in this case, which differs both 
quantitatively and qualitatively from that revealed at 
high temperatures (see Fig. 4). Therefore, more 
detailed theoretical and experimental studies are 
required for revealing the reasons for the discrepancy 
between the theory of homogeneous nucleation and 
experiment in gas-saturated solutions at low tempera¬ 
tures [20, 21], 
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